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Thermal scaling and relaxation of the interface width in an electrophoretic deposition of polymer 
chains is examined by a three-dimensional Monte Carlo simulation on a discrete lattice. Variation 
of the equilibrium interface width W r with the temperature T shows deroughening W r oc T~ s , 
with 8 ~ 1/4, at low temperatures and roughening W r oc T", with v ~ 0.4 at high temperatures. 
The roughening-deroughening transition temperature T t increases with longer chain lengths and is 
reduced by using the slower segmental dynamics. 
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I. INTRODUCTION 

Polymer chains and particulates of all shapes and sizes 
form deposits when they are driven towards a substrate. 
A continuous release and deposition of driven polymer 
chains may lead to a polymeric material with a well 
defined density at the substrate, a growing bulk, and 
an evolving interfaced The characteristics of polymer 
conformation and density at the substrate, bulk, and 
the interface depend on factors such as the nature of 
the segmental dynamics of the polymer chains, polymer- 
substrate and polymer-polymer interactions and their 
entanglement, molecular weight, temperature, type of 
medium, etc. Similar issues emerge if the polymer chains 
are driven by an electric field as in gel electrophoresis 2,3,4 
and deposit on the pore boundary of the gel. Another 
example would be a surface coating resulting from poly- 
mer chains deposited under a pressure gradient such as 
in spin-coating. The applications of such electrophoretic 
(driven) deposition and fabrication processes are enor- 
mous with a vast list of challenging questions in a rel- 
atively large parameter spaced. Using a discrete-lattice 
computer simulation model, we would like to focus here 
on the basic issue of how the interface roughness depends 
on temperature. 

We know that, as the polymer chains deposit on an 
impenetrable substrate, the polymer density increases 
and the interface develops&i&S. The fluctuation in the 
height (measured from the substrate) of the polymer 
coating along the longitudinal (i.e., growth) direction is 
the interface width W and is a measure of the surface 
roughness. The interface width is found to grow in time 
with power laws (W = At^ 1 , Bt^ 2 ) , where /3i is the ini- 
tial growth exponent followed by ( $ 2 )L2 i i2iii i i2ii£ before 
saturating to a steady-state (constant) value W s in the 
asymptotic regime as the polymer chains continue to de- 
posit. For the interface growth of particles (this could 



be viewed as the extreme limit of reducing the chain 
length) depositionifliiiiiS, the saturation of the interface 
occurs as the height-height correlation length exceeds 
the substrate length. For the electrophoretic deposition 
of polymer, such a relation between the height-height 
correlation length and substrate length is less cleari&. 
The problems become much more complex due to addi- 
tional relevent lengths such as radius of gyration, com- 
plex density-density correlations (intra- and inter-chain) , 
and their dependence on temperature, field, etc. In fact, 
we have observed interesting roughening and deroughen- 
ing of the interface as various influences (driving mecha- 
nisms such as temperature and fieldi 3 -) compete and co- 
operate. 
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Recently it was foundi 3 - that the steady-state inter- 
face width W s may reduce to a lower value (W r ) if the 
system is relaxed by continuing segmental movements of 
existing chains without adding new chains into the sys- 
tem. We have noted contrast in the characteristics of 
the relaxed interface width W r including its dependence 
on the molecular weight from that of the steady-state 
width Wg. In our previous study 1 ^, we examined the 
effect of the driving field on the interface width. De- 
pendence of the interface width and morphology on the 
temperature for particle deposition has been recently ob- 
served experimentallyi 4 * 1 ^* 1 ^ 1 ^ as well as with Monte 
Carlo simulations 1 £ii£i 1 £. In fact, the variation of the in- 
terface width with temperature is shown, under certain 
conditions, to be nonmonotonic 15,18 . Polymer surfaces 
are also affected by the deposition temperatureSiiiiS. In 
this article, we address the dependence of the interface 
width on the temperature using a variety of chain lengths, 
sample sizes and segmental dynamics after relaxing the 
polymer chains by "turning off" the source of new chains. 
In the following, we will briefly introduce the simulation 
model and follow with a discussion of the results and 
conclusions. 
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II. METHOD 

We use our previous computer simulation modcl^liSft 
on a discrete cubic lattice of size 100 xLxL with L = 4- 
40. This is a coarse-grained model with the charged 
polymer (polyelectrolyte) chains represented by a num- 
ber of contiguous occupied lattice sites, each assigned 
a unit charge. The chains are initially placed near one 
end of the sample (near x — 1) by a self- avoiding walk 
and are then driven toward the substrate/wall at the 
opposite end (x = 100) by a field E. The length of a 
chain L c is the number of lattice sites it occupies. For 
the simulations reported here, L c = 10-200. In addi- 
tion to excluded volume we consider nearest-neighbor 
and polymer-polymer repulsive interactions as well as 
polymer- wall attractive interactions. Two different com- 
binations of segmental chain dynamics 21 , kink-jump and 
crankshaft (KC) (slower)^ and kink-jump, crankshaft, 
and reptation (KCR) (faster) are used to move the chains 
according to the Metropolis algorithm. The energy 
change for a chain segment due to the field is 

AU = EAx, 




FIG. 1: Density profile using a combination of kink-jump 
and crankshaft (KC) segmental dynamics. The sample size 
is 100 x 40 x 40 with L c = 40 and E = 0.5 at various tem- 
peratures (in arbitrary units) for 10-20 independent samples. 
Density is the fraction of occupied lattice sites and x has units 
of the lattice constant. 



where Ax is the displacement parallel to the field. The 
repulsive inter- and intra-chain (Coulombic) interaction 
and attractive polymer-substrate interaction is repre- 
sented by the energy potential 

u = jY /PiPj , (i) 

ij 

where J = 1, pi = 1 for every occupied lattice site i, 
Pi = for all unoccupied sites, and p L = —1 for ev- 
ery site along the substrate (x — 101). The summation 
in eq ^ is over all nearest-neighbor sites in the lattice. 
An open boundary condition is used along the longitu- 
dinal (x) direction and a periodic boundary condition 
along the transverse (y,z) directions. An attempt to 
move each node once is defined as one Monte Carlo Step 
(MCS) . One polymer chain is released into the sample ev- 
ery 32L C /L 2 MCS for a constant deposition rate of 0.0005 
monomer units per MCS per substrate lattice site. After 
75,000 MCS, we allow the sample to relax by stopping 
the release of new chains but continuing to drive the re- 
maining chains with the electric field until equilibrium is 
reached. The simulations are performed for many inde- 
pendent samples to calculate the average values of the 
physical quantities, i.e., density profile, interface width, 
radius of gyration, etc. Further, we have used different 
sample sizes to check for the finite size effects. Most of 
the data presented here were generated on a 100 x 40 x 40 
lattice with 10-20 independent simulations, and we have 
not noticed severe finite size effects on the qualitative 
nature of the physical quantities, except for low temper- 
atures (T < 0.3). The unit for all length quantities is 
the lattice constant, and arbitrary units are used for the 
temperature and electric field. 



III. RESULTS 

Figure^shows the equilibrium density profile with the 
KC dynamics at different temperatures with E = 0.5 and 
L c = 40. Note that the shape of the density profile de- 
pends on the temperature. It is worth pointing out that 
we have also carried out our simulations at lower temper- 
ature (i.e., T < 0.1) with KC dynamics, but encountered 
technical problems, namely clogging, before deposition 
and a long relaxation time. On the other hand, the inter- 
face width (along with all other physical quantities such 
as density profile, radius of gyration, etc.) has reached 
equilibrium at all temperatures we have studied with the 
fastest (KCR) dynamics. The shape of the density pro- 
files using KCR dynamics is similar to those in Figure ^ 

The interface width is defined by 

W 2 _ Y,i(hi-hi) 2 
L 2 

where hi is the distance in units of the lattice constant 
from the substrate lattice site i to the deposited poly- 
mer farthest from the substrate (in the same y-z column) 
and h{ is the average distance from the substrate. Vari- 
ation of the interface width with time (in Monte Carlo 
Steps) for KCR dynamics is presented in Figure 2. We see 
that the initial power-law growth of the interface width 
(W oc t@) is followed by a steady-state value after 10 4 
MCS. After 7.5 x 10 4 time steps, the steady-state inter- 
face width W s is allowed to relax to its equilibrium value 
W r (by ceasing to release new polymer chains into the 
system) as in Figure 3. Results are qualitatively similar 
for the KC segmental dynamics. Note that the relax- 
ation to equilibrium occurs rapidly except at high tern- 
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FIG. 2: Interface width W (in units of the lattice constant) 
versus time steps t before surface relaxation on a log-log scale 
using KCR segmental dynamics for selected temperatures (in 
arbitrary units) and 10-20 independent samples. Other pa- 
rameters are the same as in Figure 1. 
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FIG. 3: Interface width W (in units of the lattice constant) 
versus time steps t at later times showing surface relaxation. 
The simulations used KCR segmental dynamics for different 
temperatures (in arbitrary units) and 10-20 independent sam- 
ples. Other parameters are the same as in Figure 1. 



FIG. 4: Equilibrium (relaxed) interface width W r (in units of 
the lattice constant) versus T using KC, and KCR (squares) 
segmental dynamics with E = 0.5. Different chain lengths 
are used with the KC dynamics. Other parameters are the 
same as in Figure 1. The solid line shows the power- law decay 
(cq[5J with v = —0.25 The dotted, dashed, and dotted-dashed 
lines show the power-law growths (eq0 with S =0.37, 0.38, 
and 0.25, respectively. Temperature T and E are in arbitrary 
units. 




FIG. 5: Equilibrium (relaxed) interface width W T (in units 
of the lattice constant) versus T using a variety of sample 
sizes (L = 4-40) with KC (circles) segmental dynamics and 
E — 0.5. Other parameters are the same as in Figure 1. 
Temperature T and E are in arbitrary units. 
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TABLE I: Transition Temperatures T t at E = 0.5 





L 


dynamics 


T t 


40 
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KC 


0.20 


40 


10 


KC 


0.34 


40 


20 


KC 


0.36 


40 


40 


KC 


0.35 


10 


40 


KC 


0.32 


200 


40 


KC 


0.7 


40 


40 


KCR 


3.0 



peratures (T = 4.0, 8.0). Figure 4 shows the variation of 
the equilibrium interface width W r with the temperature 
on a log-log scale using KC and KCR segmental dynam- 
ics. For KC dynamics in the high temperature regime 
(T > 4.0), W r shows a power-law scaling, 
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where v ~ 0.4. At low temperatures, on the other hand, 
there is a tendency for W r to decay with temperature. 
With KCR dynamics, we see a power-law decay, 

W r oc T~ s , (3) 

where S ~ 1/4. Thus, there is a power-law scaling for the 
roughening of the equilibrium width W r at high tempera- 
ture and deroughening at low temperatures. For L c — 40 
the roughening-deroughening transition temperature for 
KC dynamics (T t w 0.35) is lower than for KCR dy- 
namics, where T t ~ 3.0. On the other hand, the reduced 
mobility of longer chains results in an increase in the tran- 
sition temperature so that Tt ~ 1.8 for L c = 200 while 
Tt « 0.32 for L c = 10. To investigate finite-size effects we 
vary the sample size L = 4-40 and plot the temperature 
dependence of the interface width in Figure [SJ For mod- 
erate temperature (T =0.4-1.0), no significant finite-size 
effects are observed. At higher temperatures, thermal 
mobility increases the lateral height-to-height correlation 
length, which limits the interface width for small sam- 
ple sizes (L — 4). The transverse (y,z) components of 
the polymer chains' radius of gyration increases at low 
temperatures^, which leads to definite finite-size effects. 
The transition temperatures are summarized in Table [I] 

It is also of interest to compare the temperature de- 
pendence of the steady-state (W s ) and equilibrium (W r ) 
interface widths. Figure|S]shows their variation with tem- 
perature for KC dynamics. We see that the nonmono- 
tonic dependence of W s on the temperature is qualita- 
tively similar to that of W r with the faster KCR dynam- 
ics. Scaling of W s with T is perhaps more complex than 
a single power-law as for W r . 

Even for simpler particle surfaces there exists a com- 
plex temperature dependence for roughness with several 
mechanismsASi, yet the thermal scaling observed in our 
polymer interface simulations is even more complex due 
to the many degrees of freedom for the chains. At low 
temperatures we can expect an increase in temperature 
to facilitate lateral (intralayer) diffusion"^ with relatively 



FIG. 6: Interface width (W a and W r , in units of the lattice 
constant) versus T (in arbitrary units) with KC segmental 
dynamics. The parameters are the same as in Figure 1. 



less constraint due to entanglement. This will allow 
the polymer chains to reach a lower energy configura- 
tion (smoother interface). At temperatures above the 
roughening-deroughening temperature, azimuthal (inter- 
layer) diffusion dominates, and constraint due to entan- 
glement becomes more pronounced creating a (rougher) 
dynamic interface. 



IV. CONCLUSION 

In summary, the relaxation of the polymer chains, den- 
sity, and interface width depend on the segmental dy- 
namics as well as the temperature. The combination of 
slow-to-fast segmental dynamics (KCR) seems to rather 
effectively equilibrate the system in our observation time. 
Simulation with the slower KC dynamics at very low tem- 
peratures, while not very efficient for reaching equilib- 
rium, provides insight into the nature of relaxation and 
the effects of clogging. However, we do not know which 
segmental dynamics is more appropriate for comparison 
with any given laboratory measurements, though we are 
able to make some important scaling predictions: the 
equilibrium interface width W r shows roughening with 
the temperature, W r oc T 4 in the high temperature 
regime and deroughening W r oc T -1 / 4 at low temper- 
atures (at least with KCR segmental dynamics). The 
transition temperature T t that separates the roughen- 
ing and deroughening regimes depends on the polymer 
chain length and segmental dynamics. The steady-state 
interface width (W s ) also exhibits a nonmonotonic de- 
pendence on temperature, i.e., deroughening at low tem- 
peratures followed by roughening at high temperatures 
with both KC and KCR dynamics. 
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